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Abstract 

We consider the ID motion of an overdamped Brownian particle in a general potential in the 
low temperature limit. We derive an explicit expression for the probability distribution for the 
heat transferred to the particle. We find that the local minima in the potential yield divergent side 
bands in the heat distribution in addition to the divergent central peak. The position of the bands 
are determined by the potential gaps. We, moreover, determine the tails of the heat distribution. 
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There is a strong current interest in the thermodynamics and statistical mechanics of 
small fluctuating systems in contact with a heat reservoir and driven by external forces. 
The strong interest stems from the recent possibility of direct manipulation of nano systems 
and biomolecules. These techniques permit direct experimental access to the probability 
distribution functions (PDFs) for the work or for the heat exchanged with the environment 
, |9| . These techniques have also opened the way to the experimental 
verification of the fiuctuation theorems, which relate the probability of observing entro py- 



generating trajectories, with that of observing 
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consuming trajectories. [10 
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27(1 . However, from a theoretical 



point of view the time evolution of the work and heat PDFs for a Brownian particle in 
an external potential is governed by a complex partial differential equation j^, whose 
explicit solution is only available for simple potentials. It is therefore of interest to extract 
some general properties regarding the motion of a Brownian particle. 

In the present letter we consider a Brownian particle in a general static potential. We 
show that an explicit asymptotic expression for the heat PDF can be obtained in the low 
temperature - long time limit. As a starting point we consider the stochastic motion in 
ID of an overdamped Brownian particle in the general static potential U{x). Within a 
conventional Langevin equation description the motion is governed by 



dx 



dx 



(1) 



where the noise characterizing the fluctuations imparted by the heat bath at temperature T 
is correlated according to 



mm) = 2TS{t). 



(2) 



We are setting = 1 and the kinetic coefficient P = 1; note that for the motion in a viscous 
medium we have according to Stokes theorem P = l/Gvrr^i?, where R is the radius of the 
particle and t] the viscosity. 

The kinetics of a particle moving in a general potential is complex. On the one hand, 
the particle can be trapped in local minima on a time scale given by the inverse spring 
constant in a local harmonic approximation; on the other hand, the particle can also make 
a Kramers transition across potential barriers separating local minima. These transitions 
are typically dominated by the Arrhenius factor exp(— At//T), where AU is the potential 



barrier. However, in the long time limit, i.e., at times larger than the largest relaxation 
time, the particle samples the full potential profile and the stationary distribution is given 
by the Boltzmann expression 

p-PU{x) 

- (3) 

where Z{(3) is the partition function 

Z{(3) = [ dxe-^^("), (4) 



and (3 = 1/T (the inverse temperature). 

At low temperature, (3 oo, a steepest descent argument applied to the partition func- 
tion (jl]) implies that only the local minima in U contribute. The particle is trapped in 
the local minima and only rarely makes Kramers transitions to neighboring wells. Expand- 
ing the potential to quadratic order about the well at position Xi with gap Ui and second 
derivative ki (spring constant) we obtain the local contribution 

U{x)^U, + h,{x-x,y, (5) 



and for the partition function, performing the Gaussian integral, J dx exp{—ax'^) = ^Jir/a, 

m-Ei^fe-^"'. (6) 

In the present letter we wish to focus on the distribution of the heat exchange Q{t) with 
the reservoir in the long time - low temperature limit. For the Brownian motion of a particle 
in a potential the heat Q{t) sampled up to time t is a fluctuating quantity characterized by 
a time dependent probability distribution P{Q,t). In the long time limit the distribution 
P{Q,t) approaches a stationary distribution Po{Q) which we proceed to analyze. 

The heat delivered by the heat bath in the time span t is in general given by the expression 
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f dx J Q \^ dx j .j.^^f'^ dij 

where the first integral is interpreted according to the Stratonovich integration scheme and 
where we simply sample the energy change associated with the fluctuating position of the 
particle; xq is the initial position and x{t) the position at time t. The expression ([7]) can, 



of course, in principle be evaluated by the insertion of ([T]). However, this more cumbersome 
procedure can be circumvented by noting that simple quadrature yields 

Q{t) = U{x{t))-U{xo). (8) 

Note that this result only holds for a static potential. In the case of a prescribed dynamic 
potential doing work on the particle according to a given protocol, the expression ([8]) does 
not hold and one is faced with the more intractable problem of handling ([7]). 
Introducing the characteristic function the heat distribution is given by 

PiQ,t)= I ^e^?'Q(e-*fQW), (9) 

J ^TT 

where the average (■ ■ ■) in the characteristic function is both over the initial position xq and 
the final position x{t). At t = the heat distribution P((5,0) = 5{Q). After a transient 
period of order the inverse spring constants and Kramers rates the heat distribution function 
becomes stationary. The particle is at all times in thermal equilibrium characterized by the 
Boltzmann distribution ([3]). Accordingly, averaging at long times over xq and a;(oo) using 
([3]) we obtain for the characteristic function 

Finally, extending the partition function (jl]) to complex inverse temperature we obtain in a 
compact manner a general long time expression for the characteristic function 

|2 
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' Z{(3f ' ^ ' 



note that (e-^P'3(°°))p =0 = 1 ensuring the normalization condition j dQP{Q, oo) = 1. 

At low temperature the partition function predominantly samples the local minima in 
U and we obtain inserting the asymptotic expression ([6]) generalized to complex inverse 
temperature the low temperature - long time expression for the characteristic function 



|/g-«pQ(oo)^ 



^ EAhkj)-'/' exp(-/3(f/, + U,)) exp(-zp(f/, - U,)) 



+ p 



(12) 



First, we note again that (e~*^'^'^°°-*)p=o = 1 yielding the normalization of P{Q, oo), moreover, 
the phase factors exp{—ip{Ui — Uj)) can according to ([9]) be absorbed in a shift of Q. The 
interesting aspect resides in the prefactor (/9^ + p^)~^/^ which has branch points in the 
complex p plane at p = ±ip. 



By inspection of (Q and (fT2|) we note that for small Q relative to Ui — Uj the integral in ([9]) 
is logarithmically divergent for large p since ~ p~^, implying a logarithmically 

divergent contribution to P{Q, oo) = Po{Q), i.e., Po{Q) ~ — log \Q — {Ui — Uj)\. For large 
IQI, i.e., the tails of the distribution Pq{Q), we sample the small p region in ^ and f[T^ and 
closing the contour in the upper half plane (lower half plane) for Q > (Q < 0) picking up 
the branch point contribution p = i(3 {p = we obtain the dominant exponential tails 

Po(g)~exp(-/3|Q|). _ 



However, using the well-known identity 29[ | 

cos ax 

dx- 



Ko{ab), 



(13) 



(62 + x2)V2 

where -^'0(2;) is a Bessel function of the second kind, it is easy to derive an explicit expression 
for the heat distribution function. We have 



PoiQ) 



(14) 



E„^(^n,^™,)-i/^e-/'(^"+^-) 

Using the identity J dxKo{x) = it following from f|T3|) we confirm the normalization condition 
J dQPQ^Q) = 1. For small argument -^'0(2^) ~ ~ ^og{x) and we obtain for Q Ui — Uj 



PoiQ) 



-\/2^-l3{Ui+U,) 



\og\Q-{U,-U,) 



(15) 



showing that Pq{Q) exhibits a multi band structure of log divergent peaks at Q = f/j — Uj in 
agreement with our qualitative discussion. For large argument -^'0(2;) ~ (7r/2x)^/2 gxp(— a:) 
and we obtain the exponential tails Pq{Q) ~ Q"^^^ exp(— /5|(5|), including the prefactor 

In the case of the double well potential depicted in Fig. [1] with two minima, the second 
minimum with gap t/, we obtain from f|T2l) the characteristic function 
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COS pU 



(16) 



yielding the heat distribution function 

PoiQ) 



IT 
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-I'''{K^{(3\Q-U\) + Kq{(3\Q + U\)). 



(17) 
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FIG. 1: Double well potential with gap U (arbitrary units). 

In Fig. [2] we have shown the heat distribution function Po{Q) as a function of heat transfer 
Q for the parameter values ki = k2 = 1, f3 = 1, and U = 1.5. 

The interpretation of the log divergent multi band structure in the heat distribution 
'unction in the low temperature-long time limit is easy. From the work of Imparato el al. 



30l |. see also the work of van Zon et al. 



it is well-known that the heat distribution 



function for a static harmonic potential is given by the Bessel function KQ{j3\Q\) exhibiting 
a log divergence at zero heat transfer Q ~ 0; for large Q the distribution falls off according 
to the Boltzmann factor exp(— For a general potential possessing several minima 
these features seem to persist. Each minima in the potential acts like a local reservoir where 
the particle at low temperature can be trapped for a long time before making a Kramers 
transition to another well. Since the transfer between the potential well with gap Ui to the 
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FIG. 2: Heat distribution function pertaining to a double well potential with gap U. The spectrum 
shows a central peak at Q = and two side bands at Q = ±U (arbitrary units). 

well with gap Uj involves the energy difference Ui — Uj the divergent band appear at heat 
transfer Q = Ui — Uj. It also follows from (fT5l) that the contributions are weighted with the 
corresponding Boltzmann factors. 

In the case of the double well potential depicted in Fig. [1] with a well with zero gap and a 
well with gap U the discussion is particularly transparent. At low temperature the particle is 
for most of the time trapped in the zero gap well yielding the log divergent behavior for zero 
heat transfer. Occasionally the particle makes a Kramers transition to the well with gap U 
and becomes trapped yielding log divergent peaks at Q ~ ±f/. The side bands originating 
from the well with the gap are down by the Boltmann factor exp{—j3U). The sum rule 
(normalization) J dQPo^q) = 1 implies that the total integrated strength is constant. We 
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also note that for the case of a vanishing gap U = 0, i.e., for two gapless wells, we recover 



the result Po{Q) = {l3/'n')Ko{i3\Q\), independent of the spring constant [30 1. 

In order to model a double well potential in more detail we have used the fourth order 
polynomial 

Uix) = dx^ + AUx^ + - x\ (18) 

which has a gap f/, minima located at x = —{3U/2d) ± 1/2, and spring constants 
k = 2{d ± 3U). Choosing P = 1, U = 2.5 and d = 15 we have in Fig. [3] depicted the 
analytical expression (ITTl) together with a numerical solution of the Langevin equation ((10^ 
independent trajectories, with to = 0, t = 10). 

Finally, we derive a general expression for the large \Q\ behavior of the heat distribution 
function Po{Q). For large \Q\, corresponding to large heat transfer to the heat bath, we 
sample the wings of the potential. Considering a general potential U behaving like U ~ Ax^, 
for large n even, we obtain for the heat distribution function 

Po(Q) oc J da:idx2e-^(^i+^^)5(|g| - Ui + U^) (19) 

where Ui = U{xi) and U2 = U{x2)- For large \Q\ inserting U ~ Ax^ and introducing polar 
coordinates Xi = rcos0, X2 = rsin</) we have, using the delta function to eliminate r, 

PoiQ) oc [ rrfrrf0e-^l'3|F(0)^Qg| _ ^^n^^^^n ^ _ gij^n, ^^^^ ^gO) 



where F{(j)) = (1 + tan" (p) / (1 — tan" (p). For large Q the integral is dominated by the minima 
of F for = and = vr. Expanding about the minima to second order, F ^ 1 — 2(50)", 
performing the Gaussian integrals and the r-integration over the delta function in (!20l) we 
obtain the distribution function for large \Q\ 

PoiQ) oc |g|i/"-ie-^l«l. (21) 

This is a general result. Here e~^'*^' is a Boltzmann factor associated with the heat transfer 
Q whereas the prefactor IQI^^""^ is a "density of states" contribution. For n = 2 we 
obtain the previous Bessel result, Po{Q) ~ \Q\~^^'^ exp{—P\Q\), pertaining to the harmonic 
approximation, for n = 4, corresponding to the double well model potential fll8p . we have 
PoiQ) ~ IQI"^^^ exp(— /^IQI). In the inset in Fig. [3] we have depicted the heat distribution 
in a log-linear scale. The dashed line corresponds to Po{Q) oc exp(— /3|(5|). 
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FIG. 3: Analytical solution eq. ()17p (full line), and heat distribution (histogram) as obtained by 
numerical solution of the Langevin equation (10^ independent trajectories, with Iq = 0, t = 10), 
for the potential in (jlSp with d = 15, U = 2.5, (3 = 1. Inset: Heat distribution in log-linear scale. 
Dashed lines: exp(-/3|Q|). 

In this letter we have generalized the result for the heat distribution function for a har- 
monic potential to the case of a general static potential with several minima in the long time 
- low temperature limit. Our analysis shows that the gap structure of the potential wells 
give rise to a multi band structure of log divergent peaks in the heat distribution function. 
We have, moreover, derived a general result for the tails of the heat distribution function. 
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